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Abstract 

Hydrodynamic equations for a binary mixture of inelastic hard spheres are 
derived from the Boltzmann kinetic theory. A normal solution is obtained via 
the Chapman-Enskog method for states near the local homogeneous cooling 
state. The mass, heat, and momentum fluxes are determined to first order 
in the spatial gradients of the hydrodynamic fields, and the associated trans- 
port coefficients are identified. In the same way as for binary mixtures with 
elastic collisions, these coefficients are determined from a set of coupled linear 
integral equations. Practical evaluation is possible using a Sonine polynomial 
approximation, and is illustrated here by explicit calculation of the relevant 
transport coefficients: the mutual diffusion, the pressure diffusion, the thermal 
diffusion, the shear viscosity, the Dufour coefficient, the thermal conductivity, 
and the pressure energy coefficient. All these coefficients are given in terms of 
the restitution coefficients and the ratios of mass, concentration, and particle 
sizes. Interesting and new effects arise from the fact that the reference states 
for the two components have different partial temperatures, leading to addi- 
tional dependencies of the transport coefficients on the concentration. The 
results hold for arbitrary degree of inelasticity and are not limited to specific 
values of the parameters of the mixture. Applications of this theory will be 
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discussed in subsequent papers. 
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I. INTRODUCTION 



The qualitative properties of a granular gas whose dynamics is dominated by pair-wise 
collisions between the grains can be described by a Boltzmann equation, suitably modified 
to describe the more complex two particle collisions. In the simplest model the grains are 
taken to be smooth, hard spheres with inelastic collisions. In recent years, the derivation of 
hydrodynamic equations for a one component granular gas from this idealized Boltzmann 
description has been worked out in detail to Navier-Stokes order,0 with explicit expressions 
for the transport coefficients as a function of the degree of dissipation (the restitution co- 
efficient). These recent results improve upon earlier studies^ □ by providing expressions 
that are accurate even for strong dissipation. As a consequence, there are now precise pre- 
dictions from the Boltzmann equation suitable for detailed comparison with Monte Carlo 
simulation, molecular dynamics simulation, and the evolving new class of controlled experi- 
ments. This analysis for the one component system also has been extended to dense gases 
described by the Enskog equation.! Similar studies for multicomponent granular gases are 
more scarce. Existing work on multicomponent transport appears to be based on weak 
dissipation approximations. Our objective is to provide a description of hydrodynamics 
in binary granular mixtures with a comparable accuracy to that for the one component 
system, valid over the broadest parameter range including strong dissipation. The reference 
homogeneous cooling state for a binary mixture has been discussed in detail recently^ to 
provide the proper basis for analysis of transport due to spatial inhomogeneities. Those 
results are used here to describe mass, heat, and momentum transport. The expressions for 
the distribution functions, fluxes, and transport coefficients are exact to Navier-Stokes order 
(within the context of the Boltzmann equation). 

The hydrodynamic equations for a binary mixture at low density are derived from the 
coupled set of Boltzmann equations for the two species in the same manner as for a one 
component system. The solutions for the distribution functions are expanded about a local 
homogeneous cooling state that is analogous to the local equilibrium state for a gas with 
elastic collisions. The expansion is in powers of the spatial gradients of the hydrodynamic 
fields (e.g., species densities, temperature, and flow velocity) and is an extension of the 
familiar Chapman-Enskog procedure for elastic collisions. The primary technical compli- 
cation for inelastic collisions is an inherent time dependence of the reference state due to 
collisional cooling. In a one component system this occurs through the time dependence of 
the temperature defined in terms of the mean square velocity for the homogeneous cooling 
distribution. For a two component system the temperature is defined in terms of the alge- 
braic average of the mean square velocities for the two distributions. In the case of elastic 
collisions the average temperature is the same as the kinetic temperatures for each species 
in the local equilibrium state. However, a surprising result of the study in Ref. [T^ is that 
these temperatures are all different for inelastic collisions. This does not mean that there 
are additional hydrodynamic degrees of freedom since their cooling rates are the same and 
consequently, the partial kinetic temperatures still can be expressed in terms of the average 
temperature. But the relationships between these temperatures are functions of the den- 
sities for the two species and lead to a new dependence of the reference cooling states on 
these hydrodynamic variables. The consequences of this effect for the transport coefficients 
is significant, as shown below. 
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The hydro dynamic equations for a binary mixture are somewhat more complicated than 
for the one component case: six coupled equations with eight transport coefficients. The 
irreversible (dissipative) parts of the mass, heat, and momentum fluxes are calculated to 
leading order in the spatial gradients of the hydrodynamic fields. For systems with elastic 
collisions the specific set of gradients contributing to each flux is restricted by fluid symmetry, 
time reversal invariance (Onsager relations), and the form of the entropy production.0 For 
inelastic collisions only fluid symmetry applies so there is greater freedom in representing 
the fluxes and identifying associated transport coefficients. This is discussed further in 
Section [TTJ where the independent gradients are chosen to be those for the concentration 
of species one, the pressure, the temperature and the components of the flow field, with 
eight independent scalar transport coefficients. Using the Chapman- Enskog expansion the 
solutions to the Boltzmann equations are obtained to leading order in these gradients, and 
the transport coefficients are expressed in terms of the solutions to a set of coupled linear 
integral equations. 

The plan of the paper is as follows. In Sec. the coupled set of Boltzmann equations and 
the corresponding hydrodynamic equations are recalled. The Chapman-Enskog expansion 
adapted to the inelastic binary mixtures is described in Sec. [HI] to construct the distribution 
function to linear order in the gradients. This solution is used to calculate the fluxes and 
identify associated transport coefficients. A Sonine polynomial approximation is applied to 
solve the linear integral equations defining selected transport coefficients in Sec. [TV]. We get 
explicit expressions for these transport coefficients in terms of the restitution coefficients 
and the masses, concentrations, and sizes of the constituents of the mixture. Finally, the 
results are summarized and discussed in Sec. [V]. 



II. BOLTZMANN EQUATION AND CONSERVATION LAWS 



Consider a binary mixture of smooth hard spheres of masses mi and r«2, and diameters 
G\ and 02 . The inelasticity of collisions among all pairs is characterized by three independent 
constant coefficients of normal restitution an, a 22 , and «i2 — «2i> where is the restitution 
coefficient for collisions between particles of species i and j. In the low-density regime, the 
distribution functions fi(r, v; t) (i = 1,2) for the two species are determined from the set of 
nonlinear Boltzmann equations^ 

(^ + Vi-V)/i(r,vi,t) = X;^[vi|/i(t),/i(t)] • (1) 

j 

The Boltzmann collision operator [vi|/j, fj] describing the scattering of pairs of particles 



is 



Jij [vi|/i, fj) = crj J dw 2 J da Q(a ■ gi 2 )(a ■ g i2 ) 



o^ij fiir^'^^fjir^v^t) - f i (r,\ l ,t)f j (r,\ 2 ,t)\ , (2) 

where = (a* + <Jj ) /2, a is a unit vector along their line of centers, is the Heaviside 
step function, and gi2 = Vi — v 2 . The primes on the velocities denote the initial values 
{v^, v' 2 } that lead to {v 1; v 2 } following a binary (restituting) collision: 
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v i = vi - ^ji (l + a^. 1 ) (a- ■ g 12 )a-, v' 2 = v 2 + ^ (l + o^ 1 ) (a ■ g 12 )ff (3) 

where /i^- = mi/ {mi +rrij). The relevant hydro dynamic fields are the number densities rij, 
the flow velocity u, and the temperature T. They are defined in terms of moments of the 
distributions fi as 

m = J dvi/i(vi) , P u = J2j dvimiVi/i(vi) , (4) 

i 

nT = p = ^fdv 1 ^V?f i (Y 1 ), (5) 

where Vi = vi — u is the peculiar velocity, n = n\ + n 2 is the total number density, 
p = mini + m 2 n 2 is the total mass density, and p is the pressure. 

The collision operators conserve the particle number of each species and the total mo- 
mentum but the total energy is not conserved: 

/ dv 1 J ij [v 1 \f i ,f j ] = 0, (6) 
J dvirriiViJijivilfiJj] = , (7) 



]T J rfv 1 |m^ 1 2 J^[v 1 |/ i , fj) = -fnTC , (8) 

y 

where ( is identified as the "cooling rate" due to inelastic collisions among all species. From 
Eqs. @-(H), the macroscopic balance equations for the mixture can be obtained. They are 
given by 

D t m + mV ■ u + = o , (9) 

mi 

Au + ^VP = 0, (10) 

AT - - £ ^ + A (V ■ q + P : Vu) = -C T . (11) 

n m,- 3n 

In the above equations, D t = d t + u • V is the material derivative, 

}i = m t J dviVi/^Vi) (12) 
is the mass flux for species % relative to the local flow, 

P = E / dv 1 m i V 1 V 1 f i (v 1 ) (13) 
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is the total pressure tensor, and 

q = E / ^v 1 |m l V 1 2 V 1 / i (v 1 ) (14) 

i 

is the total heat flux. 

The balance equations become a closed set of hydrodynamic equations for the fields rii, u 
and T once the fluxes (|TT|)-(|T4|) and the cooling rate ( are obtained in terms of the hydrody- 
namic fields and their gradients. The resulting equations constitute the hydrodynamics for 
the mixture. Since these fluxes are explicit linear functionals of fi, a representation in terms 
of the fields results when a solution to the Boltzmann equation can be obtained as a function 
of the fields and their gradients. Such a solution is called a "normal" solution, and a prac- 
tical means to generate it for small spatial gradients is provided by the Chapman-Enskog 
met hod. Ill 



III. CHAPMAN-ENSKOG SOLUTION OF THE BOLTZMANN EQUATIONS 



The analysis of transport phenomena in fluid mixtures is considerably more complicated 
than in the case of a one component system. Not only is the number of transport coefficients 
larger but these coefficients also are functions of more parameters such as the concentra- 
tions, mass ratios, size ratios, and the three coefficients of restitution. It follows from fluid 
symmetry that the pressure tensor has the same form to first order in the gradients as for 
a one component system. As noted in the Introduction, there is more flexibility in the 
representation of the heat and mass fluxes. Even in the case of elastic collisions, several 
different (but equivalent) choices of hdyrodynamic fields are used so some care is required 
in comparing transport coefficients in the different representations. The choice here is to 
use the concentration of species 1, defined in terms of the densities by x\ = n\J {rii + n 2 ), 
together with the pressure, the temperature, and the three components of the local flow 
velocity.0 The fluxes then have the forms 



/ mim 2 n\ 
ji = - I — - — I DVx 1 



p T 



J 2 



-Jl 



(15) 



q = -T 2 D"Wx l - LVp - AVT, 



P 



a/3 



pS, 



a/3 



V ( + V a Up - -5 a pV • U ) - KS a(3 V ■ U. 



(16) 
(17) 



The transport coefficients in these equations are 

/ D \ 
D' 
D 



'p 
D" 
A 
L 
V 

V « J 



( diffusion coefficient \ 
thermal diffusion coefficient 
pressure diffusion coefficient 
Dufour coefficient 
thermal conductivity 
pressure energy coefficient 
shear viscosity 
\ bulk viscosity / 



(18) 
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For systems with elastic collisions, the thermal conductivity in a mixture is generally mea- 
sured at ji = constant, based on Onsager's relations between coefficients in j x and q. This 
is no longer an experimentally useful choice here. If in addition it is required that p = 
constant, then Vxi can be eliminated to give q in terms of jj and VT. The coefficient of 
VT in this representation is then the thermal conductivity. The bulk viscosity vanishes 
at low density, as shown below, just as in the case of elastic collisions. The objective here 
is to apply the Chapman-Enskog method for a solution to the Boltzmann equation to first 
order in the gradients, confirm the forms (|15|)-(|T7j), and determine a means to calculate the 
transport coefficients as functions of the parameters of the system. 

The Chapman-Enskog method assumes the existence of a "normal" solution in which 
all space and time dependence of the distribution function occurs through a functional 
dependence on the hydrodynamic fields 



fi(r, v x ,t) = fi [v 1 |s 1 (r,t),p(r,t),T(r,t),u(r,t)] 



(19) 



For small spatial variations, this functional dependence can be made local in space and time 
through an expansion in gradients of the fields. To generate the expansion fi is written as 
a series expansion in a formal parameter e measuring the uniformity of the system, 



(i) 



c(2) 



(20) 



where each factor of e means an implicit gradient of a hydrodynamic field. The local reference 
states fi are chosen such that they have the same moments as in Eqs. (f|) and (H), or 
equivalently, the remainder of the expansion must obey the orthogonality conditions 







I dv 1 m i v 1 /i(vi) - if } (vi 







(21) 



i 



E/dvi^I/^O-ZiVi^O, (22) 



The time derivatives of the fields are also expanded as d t = +ed^ + • • •. The coefficients 
of the time derivative expansion are identified from the balance equations (PD-flnp with a 
representation of the fluxes and the cooling rate in the macroscopic balance equations as 
a similar series through their definitions as functionals of the distributions fi. This is the 
usual Chapman-Enskog method for solving kinetic equations.El The main difference in the 
case of inelastic collisions is that the reference state has a time dependence associated with 
the cooling that is not proportional to a spatial gradient. Consequently, terms from the time 
derivative df 3 ^ are not zero. 

To zeroth order in e, the kinetic equations ([!]) become 

flf/i^E^f^- (23) 

3 

The mass and energy balance equations to this order give 

d^Xi = , T^d^T = p-'d^p = -C (0) (24) 
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where the cooling rate is determined by Eq. @ to zeroth order 

C (0) = ~E / d^-m.vlJ^fPjf] • (25) 

This homogeneous state has been studied recently^ and is discussed in more detail in the 
next Section. The time derivative in Eq. ( p3l) can be represented more usefully as 

d (o) f m = _ C ( 0) {t9t + pdp) jf) = i C (o )Vui . ( Vi/ ( 0) ) ; (26) 

where = d/dvi. The second equality follows from dimensional analysis which requires 
that the temperature dependence of fjy must be of the form 

ft^x^v^iVjvo) , (27) 

where v^it) = 2T(t)(m\ + UI2) / (mim®) is a thermal velocity defined in terms of the tem- 
perature T(t) of the mixture. The dependence on the magnitude of Vi follows from the 
requirement that to zeroth order in gradients f\ must be isotropic with respect to the pecu- 
liar velocity. The Boltzmann equations at this order can be written finally as 

k (0, V^ • (Vx/W) = E (28) 



Since the distribution functions are isotropic, it follows from Eqs. ( |T2"D and (0), that the 
zeroth order mass and heat fluxes vanish while, for the same reason, the momentum flux is 
diagonal with a coefficient that is just the sum of the partial pressures, i.e., 

i (0) = 0, q<°>=0, P^=p5 aP . (29) 

To first order in the gradients, the equation for /f is 

+ A) if + MS® = - @» + Vl • V) fP , (30) 
where it is understood that i 7^ j and the linear operators A and M. { are 

A/i (1) = - (Uft\fh + Uf^\ f? ] \ + Mf?\ ff ] \) , (3i) 

MifP = -Jij{fi°\fP]- (32) 
The action of the time derivatives <9 t on the hydrodynamic fields is 

A (1) ^i = 0, (33) 

D [% = -|V • u, (34) 
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D«T = -—V ■ u, 



2T. 



(35) 



A (1) u 



P~ l Vp, 



(36) 



where 



- (.// rU'V and use has been made of the results jj ^ = q^ -* = = 0. The 
last equality follows from the fact that the cooling rate is a scalar, and corrections to first 
order in the gradients can arise only from the divergence of the vector field. However, as 
is demonstrated below, there is no contribution to the distribution function proportional to 
this divergence. We note that this is special to the low density Boltzmann equation and 
such terms do occur at higher densities. Evaluating the right side of Eq. (|30|) gives 



(0) 



d 



n 



(0) 



Vi ■ Van - 



+ 
+ P", 



, 9x ± / p,T 

f + ~ ■ (v,/f » 



■ V hip 



2 9V! 







(0) 



1/3 



3^ Vl ' W x 



Vi ■ VlnT 



(37) 



The equation for is now 

(<9< 0) + d) + MifV = A, • Van + Bj • Vp + Ci ■ VT + D i>aP V a up. 



(38) 



The coefficients of the field gradients on the right side are functions of Vi and the hydro- 
dynamic fields. They are given by 



'A 



(o) 



(39) 



Bi(V! 



1 

V 



nT ( d 



P Wi 



f 



(o) 



(40) 



Q(v x ) = - 



i s 



2 9V! 



Vi/, 



(0) 



Vi, 



(41) 



D 







i,a(3 



= Vla^fi 



(0) 



1/? 



3^ Vl '9vT /l • 



(42) 



Note that the trace of D^ a p vanishes, confirming that the distribution function has not 
contribution from the divergence of the flow field. The solutions to Eqs. (|38|) are of the form 



fP = Ai ■ Van + Bi ■ Vp + d ■ VT + V i>aP V a up . 



(43) 
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The coefficients Ai,Bi,Ci } and T>i i<x p are functions of the peculiar velocity Vi and the hy- 
drodynamic fields. The cooling rate depends on space through its dependence on X\, p, 
and T. The time derivative acting on these quantities can be evaluated by the replace- 



ment dl — > — C*- -* (Tdr +pd p ). In addition, there are contributions from d\ acting on the 
temperature and pressure gradients given by 



dr'VT = -V (TC (0) ) = -C (0) VT - TVC 



*(o) 



2T 



VT-T 



, ^1 / P ,T P 



(44) 



a^v? = -v (k (0) ) = -c (0) vp - P v( 



(0) 



-2C (0) Vp - p 



( dC (o)\ C (o) 



(45) 



The integral equations for Ai, Bi, Ci, and V i al3 are identified as coefficients of the independent 
gradients in Eq. : 



-C (0) (Td T + pd p )+d 



Ai + MiAj = Ai+ (^Pj (pBi + Td) , 



(46) 



<(°) {Td T + pd p ) + d - 2C (0) ] Bi + MiBj = Bi + ^C l; 



(47) 



(Tdr+pd p )+Ci--C^ 



pC (0) 

^ + A*^- = Ci - ^B, : 



2T 



(48) 



-C (0) {Td T + pd p ) + d V^p + MiV itaP = D i a 



(49) 



The solutions to these linear integral equations are made unique by the orthogonality con- 
ditions (|2Tp and fl22|), i.e., 



( Ai \ 

Bi 
Ci 



, J ^vimiVi 



( ^ \ 

Bi 







(50) 



Ci 

V A,a/3 / 



(51) 
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With the functions (Ai, Bi,Ci, V^ a p) determined in this way, the solutions to the Boltzmann 
equations are determined by ( f43l) exactly to first order in the spatial gradients. 

Use of Eq. fl4~3| ) in the definitions (]T2|)-(|i4|) gives the expected forms (p)|)-(|T7l) for the 
fluxes. The transport coefficients associated with ji are identified as 

D = -^—Jdw 1 Y v A 1 , (52) 

Dp = ~^I dVlYl ' Bu (53) 
B'^-^f- [dv 1 Y 1 (54) 



3p 

The transport coefficients for the heat flux are 



D " = -^J2j l^VfV, ■ Ai, (55) 

L = -\ll J d vi \m l V 1 2 V 1 ■ B h (56) 

A = ~ / dv i ImVfVi ■ d. (57) 
Finally, the shear and bulk viscosities are 

V = --^Y1 J dv 1 m i V la V 1 pV i ^, (58) 



« = ~~ E / rfv l m i V l V i,oca = 0. (59) 

The bulk viscosity vanishes since the trace of T>^ a p vanishes, as follows from Eq. (4"9"). 

To summarize the results to this point, the solutions to the Boltzmann equations to first 
order in the spatial gradients are 

fi = // 0) + Ai ■ Vxi + Bi ■ Vp + d ■ VT + V^VaUp . (60) 

The solution to zeroth order, //° , is obtained from Eq. (p£9| ) while the functions Ai,Bi,Ci, 
and T>i >a p characterizing the solution to first order in the gradients are determined from the 
integral equations ([R])-(|59|). Calculating the mass, heat, and momentum fluxes from this 
solution one can identify the transport coefficients in terms of the integrals fl52|) -(|59"D. These 
fluxes, together with the macroscopic balance equations (P)-(|TT|), provide the closed set of 
Navier-Stokes order hydrodynamic equations for a granular binary mixture. All of these 
results are still exact and valid for arbitrary values of the restitution coefficients. 
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IV. SONINE POLYNOMIAL APPROXIMATION 



Accurate approximations to the solutions to the integral equations for /j and 
(>/4.j, 13i, Ci, T)i >a p) may be obtained using low order truncation of expansions in a series of 
Sonine polynomials. The determination of // to leading order in the Sonine expansion 
has been analyzed elsewhereS and only the main result is quoted here. The polynomials 
are defined with respect to a Gaussian weight factor. The parameters of this Gaussian are 
chosen such that the leading term in the expansion yields the exact moments of the entire 
distribution with respect to 1, Vi, and vf. The latter defines the kinetic temperatures for 
each species 

\n lTl = J dv^VffW. (61) 

For elastic collisions the temperatures Tj are the same as the global temperature T defined 
in Eq. (|5|). The condition that is normal implies that Tj is a function of rii and T, or 
equivalently 

§ = 7i(zi), (62) 
where 7« depends on the hydrodynamic state through the concentration. The determination 



of this functional dependence is worked out in Ref. [10]. The Sonine polynomials used here 
are defined with respect to Maxwellians characterized by the temperatures Tj rather than 
those characterized by T for elastic collisions. With this choice the leading deviation from 
the Maxwellians is a polynomial of degree 4 

/f = n^ 3 ^ (V*) , (63) 




1 + f (<%V? - 5W 2 + H) 



(64) 



where V* = Vi/vq, and 0j = (/%i7i) -1 - If polynomials defined in terms of T are used the 
leading correction is a polynomial of degree 2 in v *, proportional to T — Tj. The choice of 
polynomials defined in terms of T effectively resums an infinite set of terms in this second 
type of expansion. Since 7$ is a function of the concentration, a significant new contribution 
to the parameters of the integral equations for the transport coefficients occurs through 
the additional concentration dependence associated with temperature difference of the two 
species. The coefficients q in Eq. ( |64]) are determined by substitution of Eq. (|63| ) into Eq. 
(PH) and retaining all terms linear in Cj for the leading Sonine polynomial approximation. 
The reader is referred to Ref. |H] for further details. 

In the remainder of this section, we will compute the mass flux, the pressure tensor, and 
the heat flux in the leading Sonine polynomial approximations. Let us consider each one 
separately. 



12 



A. Mass flux 



Here, the lowest order Sonine polynomial approximations for A%, Bi, and Ci are obtained 
and applied to the calculation of the transport coefficients in the mass fluxes, D, D p , and 
D' . The leading Sonine approximations (lowest degree polynomial) of the quantities Ai, Bi, 
and d are 



(65) 



The coefficients {a^i, b itl , c^i} are related to the transport coefficients in this approximation 
through (H)-(|D ' 







( a i,l \ 


B> - 




hi 









hi 



c i,i 



1 

£7 
1 

1 

^7 



«2,1 



62.I 



C 2 ,l 



mim 2 n 
P 



D. 



miTi 



D 



p- 



Tn x T x 



-D'. 



(66) 
(67) 
(68) 



Here, 7 = Ti/T 2 and 5 = ni/n 2 . The coefficients {01,1; 61,1; ci,i} are determined by sub- 
stitution of Eq. (|65f) into the integral equations (ff6|)-(|48|). The details are carried out in 
Appendix |A| with the results 



a>i,i 



'A 

dxi 



InniT] 



J P ,T V dx i ) P ; 



(pb hl + Tcia) 



(69) 



61.1 



l - 1 

p V 



minT\ 



2i/ 



(70) 



c i,i 



K (o) 



2Tz/ 



where the collision frequency u is given by 
mi 



3niT 

This integral is evaluated in Appendix [B] with the result 
v (l + ai 2 N 



x 



1 - 



-V>2i7i + M1272 

1 pi (m 2 7i) 2 ci + p 2 (mi7 2 ) 2 c 2 
16p (m 2 7i + mi7 2 ) 2 



(71) 



7 J dv l V 1 • (jiafvxlA^Vi./f] -^NA^mV*]) . (72) 



(73) 
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where v e is the corresponding result for the elastic case 



and pi = miUi is the mass density of species i. These results, together with 02,1 = —£701,1, 
62,1 = —£761,1, and c 2 ,i = —£701,1 completely determine the distribution functions to first 
order in the Sonine polynomial expansion. 

The transport coefficients are identified from (|69|)-(|7TD as 



m\min \ 2 



(75) 



D> = -^—D p . (77) 

Since ji = — j2 and Vxi = — Vx2, it is expected that D should be symmetric with respect 
to interchange of particles 1 and 2 while D p and D' should be antisymmetric. This can be 
verified by noting that niTi + n 2 T 2 = nT. The expression for v agrees with the known result 
for elastic collisions.^! For the case of mechanically equivalent particles [m\ = m 2 ,au = 
«22 = «i2 = ot,on = <7 22 ) in the dilute concentration limit (p 2 — ► 0), the expression for the 
diffusion coefficient D coincides with the one recently derived in the self-diffusion problem.lll 



B. Pressure tensor 

The leading Sonine approximation for the function T>i^ is 

T>i ;a f3 — > fl,Mdi,lRi,ap(Vl), (78) 

where 

Ri, a p(Vi) = rrn (Vi a V w - ^ 2 £ a/3 ) , (79) 

and 

1 1 r 

di,i = — — ^ / dwiR^ a pDi :0l p. (80) 

The coefficients dj,i are related to the shear viscosity 77 in this approximation through Eq. 
(H) as 

V = -T 2 (ni7idi,i + ^27^2,1) • (81) 
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The integral equations for the coefficients d^\ are decoupled from the remaining transport 
coefficients. The two coefficients d^i are obtained by multiplying Eqs. (|49|) with Ri^p and 
integrating over the velocity to get the coupled set of equations 



-i i - |C (0) ria \ ( di,! 

T21 r 22 -|c (0) ; 



Tf 1 



The frequencies 7$,- are given in terms of the collision operator by 



1 1 



J d\iR i a pCi (fi : AiRi, 



a/3) , 



(82) 



(83) 



1 1 



J dViRi ;a pMi (fj, M Rj,af3) , ^ 3- 



(84) 



The evaluation of these collision integrals is given in Appendix |C|. The solution with the 
matrix elements known is elementary so that r\ can be calculated directly from Eq. (|81|). 



C. Heat flux 

The computation of the heat flux requires going to the second Sonine approximation. In 
this case, the quantities Ai, Bi and Cj are taken to be 



■Ai(Vi) - h. 



DV l + a 1>2 S 1 (V 1 ) 



A 2 (V 1 ) - h 



M 



mim 2 n 
pn 2 T 2 



DV 1 + a 2 , 2 S 2 (V 1 ) 
(85) 



1,M 



pniT] 



■A,Vi + 6i, 2 Si(Vi) 



B 2 (Vi) -> / 2 . 



P 



pn 2 Tr, 



D p V 1 + b 2t2 S 2 (V 1 ) 



(86) 



Ci(V x ) - f\ 



where 



M 



P 



Tn{T] 



D'V 1 + c lt2 S 1 {Y 1 ) 



C 2 (Vx) -> / ; 



2,M 



S<(Vi) = (Im^-^JVx. 



Tn 2 T 2 

(87) 



(88) 



In these equations, it is understood that D, D p and D' are given by Eqs. (fT5D-([T7|), respec- 
tively. The coefficients a it2 , b^ 2 and q i2 are defined as 



h2 I = 15^# / rfVlSi(Vl) ' 



V Ci 



(89) 
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Consequently, the transport coefficients appearing in the heat flux, D", L and A are given 
by 



D h 



5 /niT? n 2 7f \ . hnm 1 m 2 ( 72 , n 

o T ai,2 + 02,2 + ^r- )D, 

2 V mi m2 I 1 pi Vmi m2, 



(90) 



2 V mi m 2 ' 



5 p / 71 _ 72_ 
2 n \mi m 2 



D, 



v 



(91) 



a = - V f^i Cl2 + ^i C22 l + 5 - P (?L - ) jr. 

2 \ mi ' m 2 ' / 2 Vmi m 2 , 



(92) 



The computation of the coefficients D", L, and A is also carried out in Appendix |A|. By 
using matrix notation, the coupled set of six equations for the unknowns 



{oi )2 ; a 2j2 ; &i )2 ; & 2 , 2 ; Ci )2 ; c 2i2 } 

can be written as 

Here, X a i is the column matrix defined by the set (|93|) and A^/ is the matrix 



A 



( "11 - |c (0) 


^12 


P ( dx 1 ) p T 







\ 


^21 


^22 - |C (0) 





P ( fei ) p T 





t ( 9 5 (0> ) 








^11 - |C {0) 




-TC (0) /P 











v%\ 


V22 ~ |C (0) 





-TC (0) /P 








K (0) /2T 





^11 - C (0) 


^12 


V 








K (0) /2T 


V21 


^22-C (0) / 



Here, we have introduced the collision frequencies 

2 rrij 



15 n, 



(93) 



(94) 



(95) 



(96) 



v. 



2 rrij 



The column matrix Y is 



dviS. ; ■ Mi (AmS,-) , i ^ j 



V 2 
Y 5 



(97) 



(98) 
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where 



C^mim2nD 



3^7 (/WiJ 



-57 y dviSi • Mi(/ 2 ,wV 2 ) 



2 m\m2nD 
15 pnfT x 4 



/" <fV]Si • £l(/l lM Vi 



(99) 



C ( °V^p 1 ci 2 m lP D p 
pn\Tl 2pTi 15 pn\Tf 

-S'y y dviSi • Mi(/ 2 ,jtfV 2 



y dviSi ■ £i(/i,mVi 



(100) 



where 



C (0 V^' 2 + ci 2 mipD' 
~ TnjT? " 2TT X + 15 Tnf T x 4 

-57 y dviSi • Mi(/ 2 , tf V 2 



y dviSi • £i(/i,mVi 
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j v T/ 4 f (0) 15 



(101) 



(102) 



The corresponding expressions of the elements Y^j ^4 an d ^6 can be deduced from Eqs. 
( |100| ) and ( |101| ), respectively, by interchanging 1 <->• 2 and setting D — > D, D p - 
D' — » — .D'. The evaluation of the collision integrals ( p6|) and fl97|) is given in Appendix |D[ 
The solution to Eq. fl94j) is 



Dp and 



A 



-i 



(103) 



This relation provides an explicit expression for the coefficients a^, &j,2 and in terms of 
the restitution coefficients and the parameters of the mixture. From these expressions one 
can easily get the transport coefficients D', L and A from Eqs. (|9~0"])-(|9~2j), respectively. 



V. DISCUSSION 



The primary objective of this work has been to obtain the hydrodynamic description of 
a binary mixture of granular gases from an underlying kinetic theory. The derivation of the 
hydrodynamic equations consists of two steps. First, the macroscopic balance equations (0)- 
(P"T|) for number densities, total momentum and energy are obtained from the corresponding 
kinetic equation. Next, the fluxes and cooling rate in these equations are determined from 
a solution of the coupled Boltzmann equations given in terms of the hydrodynamic fields 
and their spatial gradients. The corresponding constitutive equations for the mass, heat and 
momentum fluxes are given to Navier-Stokes order by Eqs. (|T5|)— (|TT|) , and the associated 
transport coefficients are given by Eqs. (|52D-( ^) , fl55|) -(|57D, and (|58|)-(|59D, respectively. 
These results are exact within the context of the Boltzmann kinetic equation. A practical 
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evaluation of these coefficients is possible by means of a Sonine polynomial approxima- 
tion and the derivation and approximate results are not limited to weak inelasticity. An 
exploration of the full parameter space (mass ratio, diameters, concentrations, inelasticity 
parameters) is straightforward but beyond the scope of this presentation. We intend to cre- 
ate a library of programs for calculation of distribution functions and transport coefficients 
in one and two component granular systems, available to the public on our web site. 

The Chapman-Enskog method provides an expansion of the distribution function for 
weak spatial inhomogeneity. This means that the relative spatial variation of all hydrody- 
namic fields must be small over distances of order of the mean free path. This encompasses 
a wide range of interesting phenomena, but excludes strongly driven systems such as those 
under steady shear. For such states it is not profitable to go beyond the order considered 
here in the Chapman-Enskog expansion, but rather to use other methods not based on small 
gradients. Previous applications of the Chapman-Enskog method have typically introduced 
additional assumptions for convenience that are not internally consistent with constructing 
a solution to the Boltzmann equation. In most of these cases the reference state was cho- 
sen to be a Maxwellian, presumed to give accurate results at weak dissipation. This is not 
the has been demonstrated in the one component case.EI In addition, the effects of 

cooling in the lowest order time derivative of the temperature and pressure were neglected, 
again under the assumption of weak dissipation. Here, the reference state has been taken 
to be an exact solution to the uniform Boltzmann equation and consequently, there is no 
limitation on the parameters of the system (mass ratio, size ratio, concentrations, restitution 
coefficients). An interesting and important result of this exact analysis is that the temper- 
atures of each species must be different. Furthermore, the time derivatives are calculated 
at each order in the gradient expansion without restriction on the restitution coefficients. 
This assures that a consistent solution to the Boltzmann equation is constructed at each 
order in the gradients, without any restrictions on the system parameters. This consistency 
is reflected in the verification of solubility conditions for the integral equations determining 
the transport coefficients. These are the primary new features of this work. 

The evaluation of the transport coefficients for practical results introduces a new approxi- 
mation, truncation of an expansion for the solutions to the integral equations in polynomials. 
The leading order truncation is known to be accurate to approximately five percent in the 
case of elastic collisions. Exceptions are extreme mass ratios (e.g., electron-proton systems). 
Its validity for inelastic systems has been recently checked by Monte Carlo simulation for 
shear viscosity. The accuracy is found to be similar to that for elastic collisions. 

Recently, a seemingly similar analysis for granular binary mixtures has been given in Ref. 
TBI, also based on a two temperature description. However, this work is phenomenological 



with no attempt to solve the kinetic equation. Instead, it is assumed that the distribution 
function is a local Maxwellian. This is reasonably for estimating the dense gas collisional 
transfer contributions to the fluxes, but it predicts that all transport coefficients calculated 
here at low density should vanish. Clearly, the phenomenology is flawed. 

In a next paper, we will evaluate the expressions for the different transport coefficients 
for a variety of mass and diameter ratios. Here, as an illlustration and to give some insight 
into the influence of dissipation on transport, we consider the pressure diffusion coefficient 
D p and the diffusion coefficient D. For the sake of concreteness, consider the case an = 
a 2 2 = oci2 = a. Figure |TJ shows the reduced pressure diffusion coefficient D p (a)/D(l) as 
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a function of the restitution coefficient for an = 022, ^1/^2 = 0.25, and for two values of 
the mass ratio (mi/m 2 = 0.5 and 4). Here, D p {l) = (n x T jv e p){\ — rrixn/p) is the pressure 
diffusion coefficient for elastic collisions. We see that the deviation from the functional form 
for elastic collisions is quite important even for moderate dissipation (say a ~ 0.9). This 
tendency becomes more significant as the mass of the defect particles is larger than that of 
the excess particles. Also shown for comparison in this Fig. [I] is the result for 1711/1712 = 4 
with T\ — T2 — 1, which would be obtained if the differences in the partial temperatures 
were neglected. Clearly, inclusion of this effect makes a significant difference over the whole 
range of dissipation shown (the actual value is 7 = Ti/T 2 ~ 1.36 for a = 0.8). At the level of 
the mass flux, the main transport coefficient is the diffusion coefficient D. In Fig. we plot 
the dependence of the ratio D(a)/D(l) on a for an = 022, ^1/^2 = 0.25, and for several 
values of the mass ratio im\jm2 = 0.5, 1, and 4). As before, D(l) = (pT/mim2^ e ) is the 
diffusion coefficient for elastic collisions. The shape of these curves is very similar to those 
presented for the coefficient D p , although the influence of dissipation on D is a bit stronger 
than the one observed in the case of D p . According to the behavior of D and D p one can 
conclude that the main effect of inelasticity in collisions is to enhance the mass transport 
with respect to the elastic collisions case. 

The application of the Chapman-Enskog procedure here follows closely recent deriva- 
tions of hydrodynamics and transport coefficients for one component granular gases.01 An 
important difference in the mixture is the need for two different temperatures in the refer- 
ence local homogeneous cooling states, leading to qualitative differences in the concentration 
dependence of the transport coefficients relative to the elastic limit. This reference state is 
discussed in more detail in Ref. [H]. Effectively, the absence of energy conservation for gran- 
ular gases leads to a failure of detailed balance between the velocity distributions in the 
homogeneous state. One consequence is non Maxwellian distributions observed already for 
one component systems. For multicomponent systems a second consequence is different co- 
variances of the distributions for different species, although the cooling rates are all equal. 
The latter property implies that the different partial temperatures can be expressed in terms 
of the global temperature (as is required for a hydrodynamic description) although the func- 
tional relationships defining such partial temperatures introduce a new dependence on the 
concentrations. This leads to additional spatial gradients at first order in the Chapman- 
Enskog expansion and consequently additional contributions to the transport coefficients. 
The effect of these terms can be significant, as illustrated in Fig. [l|. Previous work on gran- 
ular mixturesi~i is restricted by weak dissipation approximations. In addition, these studies 
do not include contributions from the partial temperature differences which are significant 
even for weak dissipation. 

One interesting question is whether the mixture hydrodynamics is more or less stable than 
that of the one component case, and if such phenomena as phase separation or segregation 
can occur. Another direction of study is the extension of the present simple hydrodynamic 
state to higher densities based on the revised Enskog kinetic equation, following recent results 
for the one component fluidJl In this case the new complexity is due to the dependence 
of the collision operator on concentration through the pair correlation functions for the 
different species. Finally, we hope that the present results stimulate the performance of 
computer simulations (molecular dynamics and/or by using the Direct Simulation Monte 
Carlo method^) to study hydrodynamics in granular binary mixtures. 
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APPENDIX A: LEADING SONINE APPROXIMATIONS 

In this Appendix the coefficients {a^-; b^; qj} in the leading Sonine approximations are 
evaluated. First, let us consider the coefficients {ai,i; &i,i; c^i} determining the mass flux. 
Substitution of Eq. ( |So] ) into the integral equations (^)-(^) gives 



(Al) 



- C (0) (Td T + pd p ) oli/i^Vi + en [£i/i,mVi - <fyMi/2,wVi 

'0£(O)\ 



Ax 



dx 1 I 



fi,M (p&i,i +Tc 1A )V 1 , 



C w (T<9 r + p9 p ) - 2C w j htfiMVi + 61,1 [£i/i,mVi - iT^i/j^VJ 

jv-(o) 

= Bi H /i,AfCi,iVi, 

V 



(A2) 



(Td T + P d p )-^) 



ci,i/i,AfVi + c X) i [£i/i )M Vi - <57A^i/ 2 ,mVi 

-r i, v 

— 1 2T~ ,M ,l 



(A3) 



Next, multiplication by miVi and integrating over the velocity leads to 

< (0) (Td T + pd p ) + u] ni T iaiA = i J dV 1 m 1 V l ■ Ai + (jj^J ^ {pb x ,x + Tc 



1,1, 



(A4) 



C (0) (Td T + pd p ) - 2C (0) + v niTifei,! = ~ / dVxmxVi ■ B x + ^_ ni TiCi,i, (A5) 



j£(0) 



P 



-C (0) (Td T + pd p ) - ic (0) + v 



1 /■ K {0) 
rtiTxCi,! = - y cTVimiVi • Ci - -^rniT a & lia . 



2T 



Here, v is a collision frequency defined by 



v 



J dVxmxYx ■ [£xfx, M Vx - S-yM^M] 
-— / rfVmnVj • (JisfvxI/LivfV^/f] -57^i2[vi|A (0) ,/ 2j mV 2 ]) . 



(A6) 



(A7) 
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The self-collision terms of arising from Jn do not occur in Eq. ( |A7| ) since these con- 
serve momentum for species 1. The velocity integrals appearing in Eqs. (|A4|) - (|A6|) can be 
performed using Eqs. (|39|)-(Pl]) 



i J dVmnVi • Ax = - (^- ni T^j 



(A8) 



P,T 



1 f a\t \r -r niTl d rnmT\ 



(A9) 



^ J rfVxmxVi • d = 0. 



(A10) 



From dimensional analysis niT^a^i ~ T 1 / 2 , niTi&^i ~ T 1 / 2 /^, and ni^c^i ~ T x / 2 so the 
temperature derivatives can be performed in Eqs. (A"4)- ([A6|) and the result is 



(All) 



mT\ TC (0) 



+ Ci,i, 

p 



(A12) 



c i,i 



pC 



(0) 



2Tz/ 



6i. 



(A13) 



The solutions are given by Eqs. (|69D-([7l|) of the text. 

The analysis for the coefficients {01,2; &i,2j c i,2} (Dufour coefficient, pressure energy co- 
efficient, and thermal conductivity) is similar to that above for the mass flux. The result 
is 



-C (0) (Td T + pd p ) + u 11 ] T 3 a lj2 + T 3 u 12 a 2>2 = ^^ff J ' A ^ 



(A14) 



C (0) (Td T +pd p 



2C (0) + z/nl T 3 b h2 + T 3 u 12 b 2<2 = I r/v, Sj ■ B l2 



15 niTf 



y4^(0) 
+ Ci, 2 , 



P 



(A15) 



-C (0) (T9 r +pa p )--C (0) + ^n 



^ 3 ci, 2 + T 3 i/ 12 c 2>2 = j^ ™ 1 ^ J rfv i s i ' c i2 
-^T 2 C^c li2 . 



(A16) 
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Here, the collision frequencies Uu and z/y are defined by Eqs. (|96|), respectively and (0) and 



>-12 



Ax- 



mim 2 n 



C (0) (Td T + pd p ) -V 1 f 1>M + 



mim 2 n 



D 



x [£i(/i^V!) - ^i(/2,mV x )] - P P + D') f 1M Vu (A17) 
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-D r , 



x [A(/i,^Vi) - ^iC/a^Vx)] - C (0) ^r P' " 2£ P ) /i,mV 1; 
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Ci - ^C (0) (Td T + pd p ) ^V,h,u + jAfrD' 
rii7i i z i nil i 

x [AC/i^VO - W^VO] - ^C (0) ^r P' " AO /i,mVi, 



(A19) 



where 7$ = Tj/T. The corresponding integral equations for 02,2, ^2,2, and 02,2 can be obtained 
from Eqs. ( |A14j )— ( p\16| ) by just making the change 1 <-> 2. Dimensional analysis requires 
that a\ t 2 ~ T -3 / 2 , 6 1)2 ~ T~ 3 ^ 2 /p, and a lj2 ~ T~ 5 / 2 . When one takes into account this 
dependence and the forms of Ay, By, and Cy, one finally arrives to Eq. (H). 



APPENDIX B: EVALUATION OF u 



The collision frequency v is defined by the collision integrals in Eq. (A"7). To simplify 
the integrals, a useful identity for an arbitrary function h(vi) is given by 



/dvi/i(vi)Jy [v%\fi, fj) =°ijj dvi J dv 2 fi{r,v 1 ,t)f j (r,v 2 ,t) 

x / d<r6(tr-g 12 )(<r-gi2) \h{v'i) - h^) 



with 



vj = v x - ^(1 + ay)^ ■ §12) 



This result applies for both % — j and i ^ j. Use of this in Eq. ( jA7|) gives 
mi 



- | dViVx • (jut/i^Vi, /f] " 5 7 Ji2[/i (0) , / 2 ,mV 2 
I^^Cl + aia)-^- / dVi / cZV 2 <? 12 [/^(VO/f (V 2 )(V X • g 12 ) 

O Till 1 J J 1 



^7/i ( ° ) (V 1 )/2,m(V 2 )(V2 



,12, 



(Bl) 



(B2) 



(B3) 



where 7 = Ti/T 2 and 5 = ri\/n 2 . Substitution of the distribution functions from Eqs. (|63| 
and (HH) gives 
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Here, \i = m 1 /m 2} V* = Vi/v 



J 612 



gia/«o, 0» 



2T{m\ + m 2 )/mim 2 . In addition, the quantities X (c2, #i, 6*2) and / 6 ) 2 ) are given by 



X(c 2 ,9 1 ,0 2 ) 



1 + 



d 15' 



■GO, 



J (*i,02) = / dVJ / dV 2 *e- ei ^ 2 e- e2y "^2(Vt • g 
The integral / (#1, #2) can be performed by the change of variables 

x = Vt-V*, y = 0iVt + 02VS, 
with the Jacobian (#1 + 6 2 ) 3 . The integral becomes 



12 J 



Use of this result in (B5) and (B3) gives the desired result 

{e 1 e 2 {e l + e 2 )y 3/2 Un 2l2 e 2 + n l7 i0i) 



3 7i72 



h + 2 )' 



1 

16 



(cini7i^i^2 + c 2 n 2 >y 2 6 2 6 2 



This leads directly to the result given by Eq. (|74|) in the text. 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



APPENDIX C: EVALUATION OF t u AND Tij 

The collision frequencies Tu and Ty appearing in the evaluation of the shear viscosity are 
defined by Eqs. ( p3| ) and (0). All of these have the form 

fdv 1 V 1 V 1 J ij [v 1 \f i ,f j ] = o? j J dvi / dva/iCvO/j-Cva) 

x | d<x6(5- ■ gl2 )(5- ■ g 12 ) [V'/V/ - V1V1] , (Cl) 

where the identity ( |B1| ) has been used. This result applies for both i = j and i ^ j. Using 
the scattering rule (|B"2"|), the last term on the right hand side can be explicitly computed as 



(C2) 



yy» _ Vl Vi = -^(1 + a l3 ){a ■ g 12 ) [Gijcr + crGij 

+/iii(gi2cr + crg 12 ) - Hji{l + ay)(<7 • gi 2 )crcr] 
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Here, Gy = /i^Vi + fijiV 2 . Substitution of (|C2|) into Eq. ( |C1|) allows the angular integral 
to be performed with the result 



C Tfl ■ 

J d& 6(5- ■ g 12 )(<x • g 12 ) [V'{V'{ - Y X V X \ = - -f^Hjiil + odj) [gi 2 (G ljgl2 + g l2 Gi 



+^(3 - aii)^i2gi2gi2 - ^(1 + «ii)5 3 l 



(C3) 



Notice that the last term in ( |C3j ) vanishes when is contracted with the traceless tensor 
Rj. Now, the different collision integrals can be easily calculated by the same method as 
described in Appendix O. After a lengthy calculation, one gets 
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l_;(l_ au )a 
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In the case of mechanically equivalent particles (mi = m 2 = m, <Jij = a, a,{ 



(C6) 

a, Ci = c), 

Eqs. (p4])-(|C6|) reduce to those previously calculated in the single gas case in the determi- 
nation of the shear viscosity.0 

This completely determines 7"n and r 12 . The corresponding expressions for r 22 and r 21 
can be inferred from Eqs. ( |C4j) , (|C5|) , and ( p6|) by interchanging 1 «-> 2. 



APPENDIX D: EVALUATION OF Y u i/« AND i/, 



The collision frequencies i% z/y and Yi (i = 1, ... 6) that determine the coefficients aj i2 , 
6i )2 and Cj )2 for the heat flux are defined by the collision integrals fl9~6j), (0), (P8|), (0), ( |100| ), 
and ( 101 ). All of these have the form 
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J dv 1 S i (Vi)J ij [v 1 \f i ,f j ]=(Ti :j J dvi J dv 2 /i(vi)/j(v 2 ) 

x |d£e(<x-g 12 )(£-g 12 ) [Si(Vi') - Si(Vi)] . (Dl) 



Using the scattering rule 
puted as 



the last term on the right hand side can be explicitly com- 
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Substitution of flD2|) into (Dl) allows the angular integral to be performed with the result 
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Now the different collision integrals can be evaluated. 



1. Evaluation of K- 



The coefficients ^1,3,5 are obtained from Eqs. (P9|)- (|1Q1D . The collision integrals appearing 
in these expressions can be evaluated directly by using identical mathematical steps as before. 
After some algebra, the result is 
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lo L L 

+/4(25 - 8a 12 + 3a 2 2 ) - 2// 12 (10 + fi 21 (7a 12 - 29))" 

+^ [5 + 54/i 2 2 + 15// 2 i(ai2 - 1) + 6/^(4 - 3«i 2 + 2a 2 2 ) 

-6^12(5 + ^2i(7ai 2 - 11))] - #i# 2 [-144/4 + 2/i 12 (40 + // 2 i(49ai 2 - 95)) 

+/i 21 (45 - 35ai 2 - 2/i 21 (35 - 25a 12 + 12a 2 2 )] }} , (D4) 
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Di 2 = I dv 1 S 1 -J 12 [f 1>M y 1 ji 0) ] 



--minin 2 \/7T^2i(l + "12)0^0 (#i + #2) 1/2 



h02 



-3/2 



x {5(2/? 12 + 9 2 ) + // 21 (^ + 2 ) [5(1 - a 12 ) - 2(7a 12 - U)^ 1 
+ lSp 2 12 9^ + 2/4 (2a\ 2 - 3a 12 + 4) O^fa + 9 2 ) 2 - ^9 2 9-\9 l 4 



+ 



C'2 



, ± + 9 2 )- 2 9 1 {36 2 lf i 21 (l + a 12 ) [4/i 21 (l + a 12 ) - 5] 
lb L 

+9 2 2 [-15 + 54/i 2 2 - 20/i 2 i(3 + a 12 ) + 2/4(40 + 19ai 2 + 6a 2 2 ) 
+2 A t 12 (-20 + /i 21 (61 + 7a 12 ))} + 0^ [2// 12 (-5 + 7/i 21 (l + a 12 )) 
+/i 21 (-5(9 + 7a 12 ) + // 21 (38 + 62a 12 + 24a 2 2 )] }} , 



(D5) 



C u + V n = J dvxSi • { Jn[A (0) , /i,mV 2 ] + J u [/i lAf Vi, /f 
= -lOVvfmia^n^Ti/mi) 5 /^! + a u ) 1 - «iH 
In the above expressions we have introduced the quantity 

Pl2 = ^12^2 — ^2\d\- 



320 



ci(21a u -53) 



(D6) 



(D7) 



This completely determines the coefficients 5^. 3,5- The corresponding expressions for ^2,4,6 
can be inferred from Eqs. ( p5| )-([D7|) by interchanging 1 <->• 2. 



2. Evaluation of va and v. 



The collision frequencies uu and Uij are defined by Eqs. 



and 



. These collision 



integrals are evaluated in the same manner as those for Yj. The result is 

/ dViSi ■ Ji 2 [/i {0) , / 2 ,mS 2 ] = J dv 1 S 1 ■ J 12 [f[°\ f 2M ^V 2 N 2 - |T 2 V 2 ] = T n - |r 2 C 12 , 

(D8) 



T X2 = J dv 1 S 1 ■ J 12 [fi°\ /2,M^V2 2 V 2 ] 

= \m 1 n 1 n 2 \fKiJ, 21 {l + a 12 )a 2 2 vl (#i + d 2 y z/2 (9 1 9 2 y i/2 

o 

x {2/i 2 i^ 2 (^i + 2 ) 2 (2a 2 12 - 3a 12 + A) (80! + 50 2 ) 

-/i 2 i(^i + 6 2 ) [2(3 12 9 2 2 (89 1 + 50 2 )(7a 12 - 11) - 2# 1 fl 2 - 1 (29a 12 - 37) + 25(1 - a 12 ) 

+18f3 2 2 9 2 2 (89 1 + 50 2 ) - 2p 12 9 2 \669 1 + 250 2 ) 

+59 1 9 2 \69 1 + 110 2 ) - 5(9 1 + 9 2 )9 2 1 (69 1 + 50 2 ) 

+^{61 + 9 2 )- 2 {m\ [5 + 72^ 12 + //i 2 (-66 + /i 21 (137 - 7« 12 )) 
16 L L 
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-2/x 2 i(34 + a 12 ) + A*2i(68 - a 12 + 3a 



12/ 



+5#2 5 + 54/i 2 2 - 15/i 2 i(l - ol 12 ) + 6/i2i(4 - 3a 12 + 2a 2 2 ) - 6/^ 12 (5 + ^ 2 i(7a 12 - H) 



+20 2 2 -170 + 504/i 2 2 + // 21 (55ai2 - 17) + 2/4(151 - 62ai 2 + 39a 



12- 1 



-8/i 12 (9 + 7/221(5^12 - 13))] - 0x0| [20 - 936/4 + M2i(251 - 217ai 2 ) 
+/4(-446 + 322ai2 - 168a 2 2 ) + 2/i 12 (234 + /i 2 i(329ai 2 - 607))] } , 



(D9) 



J dvA • J 12 [/i,mS 1; / 2 (0) ] = / dviSi • J 12 [/i, A f ^t?Vi - fTiVi, /f] = W 



12 



(D10) 



■(0)1 



--mi7ii7i2V^A*2i(l + ol 12 )o\ 2 vI (6»i + ; 



\-3/2 



v-3/2 



+ 2 ) 2 (2«i 2 - 3a 12 + 4) (50i + 80 2 ) 

2/3 12 07/ 2 (50 1 + 80 2 )(7ai2 - 11) + 20 2 0r 1 (29ai2 - 37) - 25(1 - a 12 ) 



x {2^ 1 1 " 2 (0 1 

-/i 2 i(0i + 2 ) 
+ 18/? 1 2 2 0r 2 (50i + 80 2 ) + 2/?i 2 0r 1 (660i + 250 2 ) 
+50 2 1 - 1 (110 1 + 60 2 ) - 5(0i + 2 )0r 2 2 (60i + 50 2 ) 

+ ^(0i + 2 )~ 2 {l50i 3 /i 2 i(l + «i 2 )(4/i 2 i(l + a 12 ) - 5) 
lb l 

+20| [45 + 540/4 + 16^21 (au - 36) + 4/4(134 + 5a i2 + 6a 2 2 ) 
-4/x 12 (148 + // 2 i(7a 12 - 263))] + 2 2 [-30 - // 2 i(267 + 217«i 2 ) 
+14/4(17 + 29ai2 + 12a 2 2 ) + 10/xi 2 (7/i 2 i(l + a 12 - 5)) 

+0i0 2 [-315 + 270/4 - 2^21(55012 + 57) + /4(440 + 326ai 2 + 156a 2 2 ) 
+2/xi 2 (-2 + /i2i(7a 1 2 + 277))]}, 



(Dll) 



J dviSj ■ { Ju[A (0) , /i,mSi] + Jiif/i^Si, /i (0) ]} = -8 v ^n 2 a 1 2 2 miT 1 (T 1 /mi) 5 / 2 (l + a n ) 



x 



33 1 

1 H (1 - an) H ci(19 - 3an' 

16 V 7 1024 v 



(D12) 

In these expressions, C\ 2 and V 12 are given by Eqs. ( P4| ) and ( |D5|) , respectively. 

In the case of mechanically equivalent particles, the expression ( |JJ) 1 2| ) coincides with the 
one previously obtained in the context of determining the thermal conductivity in a one 
component granular gas. 0From Eqs. (jDp- flDlg) , one easily gets the expressions for v 22 and 
v 2 \ by interchanging 1 «-> 2. 
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FIGURES 



FIG. 1. Plot of the reduced pressure diffusion coefficient D p (a) / D p {l) as a function of the 
restitution coefficient a = an = Q22 = &12 for an = 022 = ci2 5 a concentration ratio n\/n2 = 0.25, 
and two different values of the mass ratio: m\/m2 = 0.5 and m\/m2 = 4. The dashed line refers 
to the case m\/m2 = 4 by assuming the equality of the partial temperatures 7 = T1/T2 = 1. 

FIG. 2. Plot of the reduced diffusion coefficient D{pt)/D(V) as a function of the restitution 
coefficient a = an = "22 = &12 for an = &22 = 0"i2> a concentration ratio n\/ri2 = 0.25, and three 
different values of the mass ratio: m\/m2 = 0.5, m\jmi = 1, and m\/m2 = 4. 
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